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Abstract 

We investigate the generalization of the spin-boson model to arbitrary spin size. 
The Born-Markov approximation is employed to derive a master equation in the 
regime of small coupling strengths to the environment. For spin one half, the master 
equation transforms into a set of Bloch equations, the solution of which is in good 
agreement with results of the spin-boson model for weak ohmic dissipation. For 
larger spins, we find a superradiance-like behavior known from the Dicke model. 
The influence of the nonresonant bosons of the dissipative environment can lead to 
the formation of a beat pattern in the dynamics of the z-component of the spin. 
The beat frequency is approximately proportional to the cutoff uj c of the spectral 
function. 



1 Introduction 



The influence of a dissipative environment on the behavior of a two-state 
system is often studied within the spin-boson model [1]. This model applies 
to many systems from various fields of physics such as a magnetic flux in 
a SQUID, electrons tunneling in chemical systems or double quantum dots, 
and two- level systems in glasses [2]. Naturally, there are other systems which 
are described by a spin greater than one half. Nuclear spins constitute one 
example. The elements gallium and arsenic used in the majority of modern 
solid state experiments have a nuclear spin of 3/2. The nuclear relaxation 
process measured in NMR experiments [3] is determined by the interaction 
with a dissipative environment, for instance a two-dimensional electron gas [4]. 
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Recent experiments have also been performed on molecular magnets. These are 
small clusters of a few atoms embedded into a crystal, which can be described 
as large spins, the most prominent examples of which (Mn 12 and Fe 8 ) are 
believed to have a total spin of 10 [5,6]. 

Large spins also describe dissipation-induced collective effects of an ensem- 
ble of identical two-level systems. The coupling to the common dissipative 
environment introduces an indirect interaction between the otherwise inde- 
pendent systems. A pseudo-spin - formally the sum of all spin one halfs - 
gives information about the degree of polarization of the ensemble. 

In this work, we employ a master equation for the description of the large 
spin derived within the Born-Markov approximation. The validity of the lat- 
ter seems to be sometimes controversial in the literature, and our aim is to 
demonstrate that a careful derivation in the exact eigenstate basis of the co- 
herent system indeed yields reliable results down to zero temperature for all 
parameter values, as long as the dissipative coupling (parameter a) is small. 



2 Model Hamiltonian 

We take the examples of section 1 as a motivation to study the generalization 
of the spin-boson model to larger spins, referred to as the large-spin model in 
the following. The Hamiltonian is given by 

H = e J Z + 2T C J X + J Z J2 l q ( a l + a - q ) + J2 UJ q a l a q- i 1 ) 

q q 

The Ji are components of a spin vector J of arbitrary size J. For J = 1/2, the 
model reduces to the spin-boson system. The bias is given by e and T c accounts 
for tunneling between adjacent eigenstates of J z . The dissipative environment 
is modeled by a bath of harmonic oscillators with creation operators a} for 
a boson in mode q and interaction strength 7 g to the spin. The collective 
character of the model is due to the spin algebra which leads to non-constant 
tunnel rates for transitions between the different states of the system. 

The generalization of the spin-boson model to a dissipative multistate system 
has been studied by several authors [7,8,9,1]. There, a particle tunnels between 
different sites with a constant tunnel rate, equivalent to a tight-binding model. 
Another possible generalization of the spin-boson model arises if also excited 
states in a double- well potential are considered [10,11]. All these systems do 
not show any collective behavior by definition. 

The Hamiltonian (1) also differs from the above cited generalizations of the 
spin-boson model with respect to the interaction term. No counter term in 
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the form of Jf is included in equation (1). Such a term does not occur if 
e. g. collective effects of an ensemble of two-level systems are described by the 
large spin. However, it is interesting to note in this context that a quadratic 
term, J|, can appear in the opposite limit, the strong-coupling regime, of the 
large-spin model. This is the case if a polaron transformation is applied to the 
Hamiltonian (1). That term has considerable consequences on the dynamics 
of the large spin [12]. 

The Hamiltonian of a large spin with dissipation is closely related to the 
Dicke-model, described by the Hamiltonian 

#Dickc = UJB Jz + JxJ2 7? ( a \ + a -q) +J2^q a \ a q- ( 2 ) 

q <? 

The Dicke Hamiltonian is frequently applied in the field of quantum optics, 
where it describes an ensemble of independent two-state atoms interacting via 
the common radiation field. This model exhibits the effect of superradiance 
as pointed out by Dicke in the original work in 1954, Ref. [13]. The decay 
of an ensemble of initially excited atoms is not exponential anymore, as one 
would expect from independent atoms. Instead, the time of the decay decreases 
inversely with the number of atoms and the maximum intensity of the emitted 
radiation increases with the square of the number of atoms. The Dicke model 
does not include direct, electrostatic tunnel coupling terms. The coupling to 
the environment is offdiagonal in contrast to the spin-boson model where 
the coupling is diagonal, i.e. to the z-component of the spin. For zero bias, 
e — 0, the large-spin Hamiltonian reduces to the Dicke model, though in a 
rotated frame of reference. Thus, the ^-component of the large-spin model 
shows superradiant behavior in an unbiased system. We will find that a similar 
statement applies for the z-component of the large spin in a biased system. 



3 Master Equation in Born-Markov Approximation 



We consider the regime of weak interactions between the large spin and the dis- 
sipative environment. The Born-Markov approximation is employed to derive 
a master equation for the density matrix of the spin. This method is pertur- 
bative in the system-reservoir coupling. In second order, the master equation 
for the density matrix p(t) of the spin reads [14] 



p{t) = - fdt' Tr Rcs 



V(t),\v(t'),p(t)®Ro 



(3) 



The interaction part of the Hamiltonian, Eq. (1), is denoted by V, the equi- 
librium density matrix of the environment by Rq, and the tilde indicates the 
interaction picture. The advantage of a master equation in this form is that 
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only operators refering to the spin system enter. The degrees of freedom of 
the reservoir are already traced out with the understanding that backaction 
effects on the bosonic bath can be neglected. 



In (3), it is assumed that the spin is brought into contact with the environ- 
ment at time t — 0. Hence, the initial density matrix of the system factorizes 
in a spin and a reservoir part. For finite times, t > 0, correlations between 
the subsystems evolve. These correlations are neglected on the right hand 
side of equation (3) in the second order Born approximation as they lead to 
effects of third or higher order in the coupling. Moreover, the Markov approx- 
imation assumes that the kernel of the integration in Eq. (3) given by terms 
like TrR es {V (t)V '(£') Rq} decays on a timescale much shorter than the typical 
timescale of the spin density matrix p(t). This is reasonable for weakly in- 
teracting systems since the dynamics of the spin in the interaction picture is 
solely caused by the coupling to the environment. Non-Markovian effects in 
the Born approximation for J = 1/2 have been discussed recently by Loss and 
DiVincenzo [15]. 

Inserting V(t) in the master equation (3) and transforming back into Schrodin- 
ger picture yields the final form of the master equation for a large spin, 



p{t) = i [p(t),eJ z + 2T C J x ] - ^ (e 2 T + 4T C 2 T c ) [j z , J z p(t) 



2T c e 
A 2 
1 



(r-r c ) j z ,J xP (t) 



2T C „ 



Jz, Jyp{t) 



f — (e 2 r* + 4T c 2 r:) [J z ,p(t)j z 



(4) 



+ 



2T c e 
A 2 



(r-r*) J z , P (t)j x 



2T, 



A s 



J Z ,p(t) Jy 



Here, the influence of the environment is expressed by the rates 

r c = lp(A)coth(f)-l/> pH (^ + ^> 

T s = I fdu p{u) coth(^) ( -J^ - i I P (A), 

2 Jo yy ' \ 2 J \uj + A ou-AJ 2 Hy h 



and r = r c ( A — > 0) , where A = y 4T 2 +e 2 is the level spacing of the unperturbed 
large spin. In the following, we focus on the case of an ohmic dissipation 
when the spectral function p(ou) is linear with an exponential cutoff, p(ou) = 
2auexp(—uj/uj c ). Then, analytic expressions can be obtained for the rates 
except for the real part of T s which remains to be calculated numerically at 
finite temperatures. 
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4 Relation to the Spin-Boson Model 



For the smallest possible value of the spin, J— 1/2, the large-spin Hamiltonian, 
Eq. (1), reduces to the spin-boson model. This enables us to compare the Born- 
Markov approximation with results for the spin-boson model in the literature. 
In particular, we employ the approximate solutions of the latter as given by 
Weiss [1]. 

It is a priori clear that the Born-Markov approximation being perturbative in 
the spin-environment coupling should only be applied in the regime of small 
couplings. Within the Born-Markov approximation, the spin is described by 
the master equation (4). A peculiarity occurs for spin one half: In that case, 
a closed set of equations for the expectation values of the spin components 
follows from the master equation, comparable to Bloch equations, 

+ ^ (Re{r} - Re{r c }) (J,) + | Im{r s }, 
(J y ) = e (J x ) - ± (e 2 Re{r} + 4T C 2 Re{r c }) (J y ) ( 6 ) 
- (2T C + ^Re{r s }) (J z ) + ^| (im{r} - im{r c }), 

(J z ) = 2T c (J y ). 

A similar set of Bloch equations was recently published by Hartmann et al. [16]. 
It is, however, not possible to apply these equations to larger spins without 
additional approximations. This was already remarked in the original works 
by Bloch [17,18]. The difficulty arises because products of spin operators can- 
not be replaced by a single spin operator for spins larger than one half. The 
corresponding equations do not form a closed set anymore. Hence, we can 
only transform the master equation into Bloch equations for spin one half. 
For larger spins we will have to come back to the master equation (4). 



The equilibrium values of the spin components follow from the Bloch equa- 
tions (6). It is readily seen that {J y ) ^ vanishes. The other two components 
become 

m g Re{e 2 r + 4T C 2 T c } Im{r-r c } + eA 3 ImlTj 

{ z)co 2A Re{e 2 V + 4T C 2 T C }(A + Re{r s } - 2e 2 A 2 Re{r-r c } ' 

(7) 

m T C A 2 Im{r a }(A + Re{rj) + e 2 T c Re{r-rj im{r-r c } 
{ x)o ° ARe{e 2 r + 4T 2 r c }(A + Re{r s }) -e 2 A 2 Re{r-r c } ' 

In the limit of zero coupling, a — > 0, we retrieve with Eq. (5) the thermody- 
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namic expressions 




oo 



= — tanh 




The dynamics of the spin components follows from the numerical integration 
of the Bloch equations (6). 

For a weak ohmic dissipation and J = 1/2, Weiss gives an approximate so- 
lution for the dynamics of the spin [1]. Two temperature regimes have to be 
distinguished: At intermediate temperatures, the solution is obtained by the 
noninteracting-blip approximation (NIBA). At low temperatures, however, the 
NIBA breaks down and the solution is derived by taking into account interblip 
correlations. In the NIBA regime, the solution is given by equations (21.132) 
and (21.134) in Ref. [1]. The solution for low temperatures follows as equa- 
tion (21.172) and (21.173) in the same reference. Minor differences in the 
notation have to be remarked: The bias e, the cutoff cu c , and the tempera- 
ture ksT are identically defined in this article and in Ref. [1]. The coupling 
strength a corresponds to K in Ref. [1]. The tunnel rate typically referred to 
as A in the spin-boson literature is expressed in our notation by 2T C (here, 
A is defined as the level spacing of the unperturbed spin). Due to a different 
sign of the tunnel term in the spin-boson Hamiltonian and the Hamiltonian 
for the large spin, Eq. (1), we have (a z ) = 2 (J z ) and (a x } = —2 (J x ). 

Figures la and lb show the time evolution of J z and J x for weak ohmic dis- 
sipation, a = 0.05, at a finite temperature ksT = 2T C . The dynamics of J y 
follows as the time derivative of J z , as can be seen from the Heisenberg equa- 
tion of motion. For that temperature, the noninteracting-blip approximation 
applies and its solution is in excellent agreement with the solution of the Bloch 
equations (6). Deviations between the two methods become visible for a larger 
coupling strength, a = 0.2, as shown in the inset of the same figure. However, 
this is not surprising as the Bloch equations are perturbative in the coupling 
strength a and thus limited to small couplings a. Figure lc and Id show the 
dynamics of an unbiased spin at zero temperature. In that regime, the NIBA 
is not valid anymore and the solution of the Bloch equations is compared 
to Eqs. (21.172) and (21.173) of Ref. [1]. Again, both solutions are in good 
agreement. For smaller couplings, e. g. a = 0.025, the different solutions for J z 
cannot be distinguished anymore (not shown). The equilibrium value (Jx)^ 
of the Bloch equations appears too large as compared to the low temperature 
solution of the spin-boson model. However, both equilibrium values approach 
the thermodynamic result, Eq. (8), in the limit of zero coupling and hence 
coincide in that limit. 

We conclude that the Born-Markov approximation correctly describes the dy- 
namics of the spin-boson model for weak ohmic dissipation at all temperatures. 
This is corroborated by recent results for the driven two-state system [16]. As 
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t/T- 1 t/T- 1 

Fig. 1. Dynamics of a spin 1/2 according to the Bloch equations (6) (solid line) and 
the approximate solutions by Weiss (dashed line); a and b: e = T c , a = 0.05 (inset: 
a = 0.2), uj c = 50T C , and k B T = 2T C (dashed line: Eqs. (21.132) and (21.134) of 
Ref. [1]); c and d: e = 0, a = 0.05, lo c = 50T c , and k B T = (dashed line: Eqs. (21.172) 
and (21.173) of Ref. [1]). 

we see no reason why the Born-Markov approximation should break down for 
larger spins, we expect that the master equation (4) gives a reliable description 
of a large spin with weak dissipation. 



5 Dicke Effect 



It was pointed out in the introduction that the large-spin Hamiltonian, Eq. (1), 
with zero bias can be mapped on the Dicke model, Eq. (2), by rotation around 
the (spin) y-axis, 

#Dkkc = e in/2 J yH e~ in/2 J y . (9) 

The effect of superradiance characteristic for the Dicke model is therefore 
visible in the x-component of the large spin. This becomes apparent if the 
initial values are chosen such that (J x ) = J- Then, the effect results in an 
accelerated decay of J x with increasing spin size J. Typically, the main interest 
lies in the z-component of the spin and this is chosen maximum as initial 
value. Consequently, the question arises if that component also decays in a 
superradiant fashion. This is indeed the case as will be shown in the following. 
We choose a large bias, e = 10T C , to ensure the decay of J z . The dynamics 
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Fig. 2. Decay of J z for different spin sizes (e = 10T C , a = 0.005, to c = 50T c , k B T = T c ). 

for different spin sizes J follows from the numerical solution of the master 
equation (4). The decay of the normalized ^-component (J z ) / J is plotted in 
figure 2 for the spin sizes J = 1/2, 2, 5, and 10. It is clearly visible that the 
time in which the spin decays decreases with increasing spin size, similar to 
the Dicke superradiance. 

Naturally, the exact dynamics of the large spin differs from the Dicke model 
as the Hamiltonians are not identical. In contrast to the Dicke model, the 
coupling to the environment is diagonal and a tunnel term exists in the un- 
perturbed Hamiltonian. Yet, despite these differences, the basic results are the 
same. The reason is that superradiance is driven by the spin algebra on which 
both models rely. This becomes obvious by considering the state dependent 
transition matrix elements [13] which are maximal for the state | J, M — 0) , 
corresponding to (J z ) = 0. From the point of view of the spin-boson model 
it appears that the spin shows collective effects, namely a superradiance like 
behavior, once it is generalized to spins larger than one half. 



6 Quantum Beats 

We shall return once again to the symmetric system, e = 0. For a spin one half, 
the Bloch equations (6) predict an exponential decay for J x and a damped os- 
cillation for J z . Already for the next higher spin, J — l, new features appear in 
the dynamics of the spin. Figure 3 shows the time evolution of the expectation 
values of J z and J x for different interaction strengths. A clear beat pattern 
is visible on top of the damped oscillations of J z . With increasing interaction 
with the environment, the beat pattern dissolves to a seemingly chaotic be- 
havior. The decay of J x is superposed with oscillations. These new features 
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Fig. 3. Time evolution of J z (black line) and J x (grey line) for different interaction 
strengths with the environment (J=l, e = 0, lu c = 50T c , and /cbT = 0). 

have different origins. The beat pattern of J z is caused by the nonresonant 
bosons of the environment. The oscillations of J x , on the other hand, are due 
to two-boson processes. The latter is best understood in a rotated frame, i. e. 
in the Dicke model. There, operator combinations J+J+ and J_J_ appear in 
the master equation. In standard approaches to superradiance [19], these are 
often disregarded by invoking a secular approximation. One can show, how- 
ever, that these terms lead to an additional oscillation the frequency of which 
is approximately given by 4T C , that is twice the level spacing of the unbiased 
spin. We numerically find that the amplitude of these oscillations is increased 
by the nonresonant bosons. 



The beat pattern in the dynamics of J z is solely caused by nonresonant bosons. 
They lead to different corrections to the eigenenergies of the spin, lifting the 
equidistance of the spectrum. We consider a spin one with zero bias at zero 
temperature, the same parameters as used in Fig. 3. The three eigenstates 
of the unperturbed spin labeled in the following by |+), |0), and |— ) have 
the eigenenergies 2T C , 0, and — 2T C , respectively. The expectation value of J z 
becomes (J z ) = cos(2T c t). The nonresonant bosons of the environment lead 
to corrections to these eigenenergies. In second order perturbation theory, we 
find 

' ' (10) 



h \±) 



1 f°° 

~2 Jo dujp ^ 



F( 2 ) 



p(2) , p,(2) 



As a consequence, the three states are not equidistant anymore. The resonant 
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bosons, oo q = 2T c , are neglected in this calculation. They lead to an imaginary 
part in the energies resulting in damping, which is not our concern at this 
stage. The next step is to determine the effect of the corrections on the dy- 
namics of J z . Neglecting the corrections to the eigenstates themselves yields 

(Jz) t — COs(LU t) COs(Ld b t) (11) 

with u and u h defined as 

uj = 2T C + l - (E\% - E\%), u b = - l - E$. (12) 

As a result, the oscillation frequency u is slightly changed with respect to the 
unperturbed case and we find indeed a beat pattern with frequency uj h . For 
an ohmic dissipation, the beat frequency becomes 



cob = ot uj c + a T c 



\ u c J \uo c 



(13) 



where Ei(x) denotes the exponential integral. The term in the brackets is of 
the order of T c /u c . Hence, for u c ^> T c , the beat frequency is given in good 
approximation by uj, = «w c . For weak interactions, the resulting dynamics is 
in excellent agreement with the numerical solution of the master equation. 
Naturally, equation (11) does not include damping since dissipative effects of 
the resonant bosons are not considered in the derivation. The occurrence of a 
beat pattern is not restricted to the parameters chosen in this example. Similar 
patterns are equally found for larger spins, J> 1, and finite temperatures. The 
pattern dissolves at high temperatures or high spins. Even for a finite bias, 
e>0, such patterns can be observed. 



7 Conclusion 



The collective character of the model becomes visible for larger spins where 
we found a superradiance-like decay of the ^-component of the spin, the time 
scale of which decreases with increasing spin size. The beat pattern in the 
dynamics of J z for intermediate spins is caused by the nonresonant bosons 
of the environment. We conjecture that the beat pattern is a characteristic 
property of an intermediate spin with weak dissipation. 

In our study of a large spin coupled to a dissipative environment, the theoreti- 
cal description applied to real spins as well as to dissipation-induced collective 
effects in pseudo-spin systems such as ensembles or arrays of two-level systems. 
We derived a master equation for arbitrary spin J within the Born-Markov 
approximation. When comparing our results for J =1/2 and ohmic dissipa- 
tion to standard spin-boson NIBA (and beyond) calculations [1], we found 
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good agreement at all temperatures and arbitrary parameter values, as long 
as the dissipative coupling is weak. We conclude that the master equation in 
Born-Markov approximation is a reliable method to describe the physics in 
the weak coupling limit. 
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